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Introduction 


Point Spread Function 

- A method of evaluating the spatial 
resolution of an imaging system. 

- A measure of the spread of a single; 
point of light. 



Modulation Transfer function 
(MTF) 

- MTF is a measure of the spatial 
frequency response. 

- MTF is often calculated from the point 
spread function (PSF). 

- System response at the Nyquist 
frequency (or 0.5 cycle/pixel) is often 
used as a figure of merit. 


MTF 



H{a> x ,co y ) = 3{PSF(x,y)} 


MTF(co x , to ) 


H(<o x , to, ) 

H( 0 , 0 ) 


2-dimensional PSF and MTF are difficult to obtain. 

Often 1 dimensional functions are used: 

— 1-D PSF is the line spread function (LSF). 

- LSF can be obtained by differentiation of the edge spread function 
(ESF). 
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Techniques 


Edge Method 

- Sub-pixel edge locations were found by Fermi function fit. 

- A least-square error line was calculated through the edge locations. 

- Modified Savitzky-Golay filtering was applied on each line. 

- The filtered profile was differentiated to obtain LSF 

- MTF calculated by applying Fourier transform to LSF. 


Test image 
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Edge Method 
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• Pulse method 

- A pulse input is given to the imaging system. 

- Output of the system is the resulting image. 

- Edge detection and mSG filtering was applied to obtain 
output profile. 

- Take Fourier transform of the input and output. 

- MTF is calculated by dividing output by input and 
normalizing DC component to unity. 


nput 



Pulse method 
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SNR Definition 



_ DN difference 

SNR = 

(STD_bright + STD_dark) / 2 


Averaged Spline with data points 



DN difference 

SNR = 

STD background 
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Target Description 



Field Plan 
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Field campaign pictures on 6-22-2005 







NASA Stennis Tarp Target 

- Radiometrically and spectrally stable target with a large 
DN difference from 3.6% and 52.1% reflectance 
panels. 

- Edge oriented to obtained sub-sampled edge profile. 

- Blue tarps oriented at same angle. 



True North 

| 6° Edge Angle 



NASA Stennis tarps 
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IKONOS Acquisitions 

• IKONOS Scene Information 


Date 

Sensor 

Targets 

Resampling or 
MTF processing 

Product Type 

8/1/2005 

IKONOS 

Stennis tarps 
Blue tarps 
Mirrors 

CC / MTFC On 

Standard Geometrically 
Corrected 

CC / MTFC Off 

Standard Geometrically 
Corrected 















Date 

Sensor 

Band 

Product 

Target 

Elevation 

Azimuth 

8/1/2005 

IKONOS 

Pan. 

CC/MTF off 

Stennis 

66.0 

118.7 


mSG filtering with raw data 




Trimmed LSF 



MTF 



Date 

Sensor 

Band 

Product 

Target 

Elevation 

Azimuth 

8/1/2005 

IKONOS 

Pan. 

CC/MTF on 

Stennis 

66.0 

118.7 



mSG filtering with raw data 
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Trimmed LSF 
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Normalized frequency 


Date 

Sensor 

Band 

Product 

Target 

Elevation 

Azimuth 

8/1/2005 

IKONOS 

Blue 

CC/MTF off 

Blue tarp 

66.0 

118.7 


CD 

X 

CL 



Pixel 


Modified S-Golay interpolation with raw data 



Pixel 


Trimmed LSF 




Continuous Fourier tranform 
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MTF from Continuous Input 








Value at Nyquist = 0.3745 
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Date 

Sensor 

Band 

Product 

Target 

Elevation 

Azimuth 

8/1/2005 

IKONOS 

Blue 

CC/MTF on 

Blue tarp 

66.0 

118.7 


Angle - 5.850 [del 31 image 



Modified S-Golay interpolation with raw data 



Trimmed LSF 




Continuous Fourier tranform MTF from Continuous Input 
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LSF 


Multi-year IKONOS Comparison 

• Consistent FWHM: 1.61+/- 0.08 

• Values at Nyquist frequency were very stable. 


- Mean = 0.11, STD = 0.01 
• No trends in PSF/MTF over 5 years. 


Date 

Sensor 

Band 

Product 

Target 

8/1/05, 7/22/02, 9/26/01, 8/5/01, 
6/22/01, 5/31/01, 3/26/00 

IKONOS 

Pan. 

CC/MTFC 

off 

Stennis / Big 
Springs 


LSF over plot for Big Spring target (2000-2001) and Setnnis tarp (2002, 2005^ 


1.6 

1.4 

1.2 

1 


3-26-00 T7: 
3-26-00 17: 

5- 31-2001 

6 - 22-2001 

8- 05-2001 

9- 26-2001 

7- 22-2002 

8- 01-2005 


:T3^FWHM 
14 + FWHM 

■ FWHM= 1 .! 

■ FWHM = 1 .! 

■ FWHM = 1 J 

■ FWHM= 1, 

■ FWHM= 1 J 

■ FWHM= l.i 


^ 1 .7597 [pixel], SNR =97.2 
1.6221 [pixel], SNR = 108.3 

$626 [ppceq,5Nf+= = 457^ 

5352 [pixel], SNR = 62.4 
6380 [pixel], SNR = 79.7 

4733 [pi)cer|,^N$l = 60.9| 

6175 [pixel], SNR = 113.6 
6226 [pixe I], SNR= 92.5 


3-26-00 17:13 
3-26-00 17:14 

5- 31-01 

6 - 22-01 

8- 05-01 

9- 26-01 
7-22-02 



MTF over plot for Big Spring target (2000-2001) and Setnnis tarp (2002, 2005) 

3-26-00 17:13 
3-26-00 17:14 



3-26-00 17:13 r atlSTyqQisf = 0.0923 
3-26-00 1 7:1 4. at Nyquist = 0.1 054 

5- 31 -01 - at Nyquist = $.1 099 

6- 22-01 - at Nyquist = 0.1 207 

8- 05-01 - at Nyquist = 0.1 034 

9- 26-01 - at Nyquist = $.1 375 

7- 22-02- at Nyquist = 0.1 1 1 9 

8- 01 -05- at Nyquist = $.1117 


IKONOS SDSU campus on 8/1/2005 



CC With MTFC 


CC Without MTFC 


Noisier & 
Sharper - 
than CC 





IKONOS Summary 

• Very consistent, high-quality sensor 

- Pan FWHM =1.61 +/- 0.08 

- Pan MTF @ Nyquist = 0. 1 1 +/- 0.01 

- Multispectral MTF @ Nyquist = 0.52 (2005 data) 

• ‘MTFC on’ processing provides increased MTF response 
with typical trade-off of increased contrast with some 
additional noise in Pan, less in Blue band. 

- Pan FWHM = 1.08 +/- 0.10 

- Pan MTF @ Nyquist = 0.48 +/- 0.08 

• No indication of sensor degradation in five years. 
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Quickbird Acquisitions 

• Quickbird scene information 


Date 

Sensor 

Targets 

Resampling or 
MTF processing 

Product Type 

8/30/2004 

Quickbird 

Blue tarps 
Mirrors 

Resampling 
Kernel 
= CC 

Standard2A 

(Radiometrically corrected) 

MTF 

Standard2A 

10/5/2004 

Quickbird 

Blue tarps 
Mirrors 

CC 

Standard2A 

MTF 

Standard2A 

6/22/2005 

QuickBird 

Vertical Stennis 
tarps 

Blue tarps 
Mirrors 

CC 

Standard2A 

MTF 

Standard2A 

10/18/2005 

QuickBird 

Stennis tarps 
Blue tarps 
Mirrors 

MTF 

Standard2A 

CC 

Standard2A 
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Quickbird 8-30-2004 



Quickbird 10-5-2004 







Date 

Sensor 

Band 

Resampling 

Target 

Elevation 

Azimuth 

6/22/2005 

Quickbird 

Pan. 

CC 

Stennis along 

88.2 

75.9 


mSG filtering with raw data 
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Trimmed LSF 
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Date 

Sensor 

Band 

Resampling 

Target 

Elevation 

Azimuth 

6/22/2005 

Quickbird 

Pan. 

MTF 

Stennis along 

88.2 

75.9 


Anale=-5.681 fdeal TeSt image 
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mSG filtering with raw data 
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Trimmed LSF 
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Along-track Stennis tarp target MTF result of QuickBird on 6/22/2005 (MTF) 



Date 

Sensor 

Band 

Resampling 

Target 

Elevation 

Azimuth 

6/22/2005 

Quickbird 

Blue 

CC 

Blue tarp cross 

88.2 

75.9 


CD 

X 

CL 


Fermi edge location 



5 10 15 

Pixel 


Modified S-Golay interpolation with raw data 



Pixel 


Trimmed LSF 



Continuous Fourier tranform 




MTF from Continuous Input 
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Date 

Sensor 

Band 

Resampling 

Target 

Elevation 

Azimuth 

6/22/2005 

Quickbird 

Blue 

MTF 

Blue tarp cross 

88.2 

75.9 


CD 

X 

Q_ 





Continuous Fourier tranform 
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Normalized LSF 


QuickBird panchromatic band Along / Cross 
Track Direction Comparison 


Date 

Sensor 

Band 

Resampling 

Target 

6/22/05, 

10/18/05 

Quickbird 

Pan. 

MTF / CC 

Stennis 


• First estimate of along-track PSF/MTF. 

• Along-track PSF/MTF not significantly different than cross- 
track. 


Stennis tarp LSF plots on 6/22/05 (CC) 10/18/05 (MTF) from QuickBird 



Stennis tarp MTF plots on 6/22/05 (CC) 10/18/05 (MTF) from QuickBird 



Normalized Pulse Response Function (PRF) 


QuickBird Blue Band Along / Cross 
Track Direction Comparison 


Date 

Sensor 

Band 

Resampling 

Target 

6/22/05, 

10/18/05 

Quickbird 

Blue 

CC only 

Blue tarps 
Cross / along-track 


• Cross-track profile exhibits narrower PRF FWHM and under shoots. 

• System MTF shape differs in the orthogonal directions. 


PRF plots for QuickBird on 6/22/05 and 10/18/05 with CC interpolated scenes 


6/22/05 Cross-T. FWHM;= 2.3216 [pixel] L SNR = 113.2 
10/18/05 Cross-T. FWHM = 2.3264 [pixel], SNR = 94.6 
6/22/05 Along-T. FWHM = 2.4979 [pixel], SNR = 72.0 
10/18/05 Along-T. FWHM= 2.4777 [pixel], SNR = 22.4 


6/22/05 Cross-T. 
10/18/05 Cross-T. 
6/22/05 Along-T. 
10/18/05 Along-T. 


MTF plots for QuickBird on 6/22/05 and 10/18/05 with CC interpolated scenes 


- 0.2 




Normalized Pulse Response Function (PRF) 


QuickBird Blue Band Along / Cross 


Track Direction Comparison 


Date 

Sensor 

Band 

Resampling 

Target 

6/22/05, 

10/18/05 

Quickbird 

Blue 

MTF only 

Blue tarps 
Cross / along-track 


• FWHM values were reduced by MTF resampling process. 

• MTF values were increased — most significantly in cross track 
direction. 

PRF plots for QuickBird on 6/22/05 and 10/18/05 with MTF resampled scenes MTF plots for QuickBird on 6/22/05 and 10/18/05 with MTF resampled scenes 



Quickbird Results Summary 2004-2005 


Band 

Pan. 

Resamo. 

CC 

Scan 

Direction 

Cross 

Along 

FWHM 

(gsd) 

MTF at Nyquist 
(cycles/gsd) 

Obs. 

1 

1 

mean 

1.89 

1.89 

siama 

mean 

0.09 

0.04 

siama 


MTF 

Cross 

1.71 


0.26 


1 



Along 

1.69 


0.10 


1 

Blue 

CC 

Cross 

2.32 

0.01 

■19 

0.07 

4 



Along 

2.49 

0.01 

B 

0.01 

2 


MTF 

Cross 

2.11 

0.05 

0.30 


4 



Along 

2.27 

0.02 

0.25 


2 


Notes: 

1 . Blue Band FWHM’s are for Pulse Response Functions. 

2. Only one Pan band observation! Essentially same PSF/MTF in 
both scan directions. 

3. MTF resampling provides noticeable contrast improvement; 
SNR is lowered, but still acceptable. 

4. Overall SNR is good to excellent. 

5. PRF is very repeatable in Blue band; MTF is noticeably less so 


28 













Ouickbird SDSU campus on 6/22/2005 



MTF resampling 


CC resampling 
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Normalized LSF 


2002/2003/2005 QuickBird Comparisons 


Date 

Sensor 

Band 

Resampling 

Target 

7/20/02, 8/25/02, 9/7/02, 
9/15/03, 10/18/05 

Quickbird 

Pan. 

CC 

Stennis tarp 


• GSD changed from 0.7 to 0.6 meters after 2003. 

• Spatial resolution performance appears consistent from 2002 to 2005 . 



MTF over plot for Stennis tarp target 



Overall Quickbird Conclusions 


• Initial along-scan PSF/MTF estimates indicate 
slightly more blur than cross-scan, as expected. 

• Good to excellent SNR. 

• MTF compensation provides noticeable contrast 
boost with normal loss of SNR in Panchromatic 
band, but minimal loss of SNR in Blue band. 

• No degradation of Panchromatic band PSF/MTF 
indicated from 2002 through 2005 . 
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